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Abstract 

We use the QCD sum rule approach to calculate the masses of the Aq and 
Sq baryons to the Aqcd/^q order within the framework of heavy quark 
effective theory. We compare the direct approach and the covariant approach 
to this problem. Two forms of currents have been adopted in our calculation 
and their effects on the results are discussed. Numerical results obtained in 
both direct and covariant approaches are presented. The splitting between 
spin 1/2 and 3/2 doublets derived from our calculation is Sg^ — Eg ~ 
0.35 lb 0.03 GeV^ which is in good agreement with the experiment. 
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I. INTRODUCTION 



Important progresses in the theoretical description of hadrons containing a heavy quark 
have been achieved with the development of the Heavy Quark Effective Theory (HQET) 
Based on the spin-flavor symmetry of QCD, exactly valid in the infinite tuq limit, 
this framework provides a systematic expansion of heavy hadron spectra and both the 
strong and weak transition amplitudes in terms of the leading contribution, plus corrections 
decreasing as powers of l/mq. HQET has been applied successfully to learn about the 
properties of mesons and baryons made of both heavy and light quarks. 

The effective Lagrangian of the HQET, up to order I/itlq, can be written as 

C,s = Ktv-DK + -^IC + -^S + 0{l/ml), (1) 

where hv{x) is the heavy quark field in effective theory. Apart from leading contri- 
bution, the Lagrangian density contains to OII/uiq) accuracy two additional opera- 
tors /C and S. K, = h^{iD^)'^h^ is the non-relativistic kinetic energy operator and 
^ ~ ^ ( °'a^/^ ^ ) ° h^'^p.vQsG^" is the chromo-magnetic interaction. Here {D^Y = 

D'^ — [v ■ -D)^, with Dfj^ = dfj_ — i gA^ the covariant derivative and (3q = 11 — ^Uf is the 
first coefficient of the f3 function. 

The matrix elements of the operators /C and S in (|l|) play a most significant role in 
many phenomenological applications such as the spectroscopy of heavy hadrons IQ and 
the description of inclusive decay rates 0. For the ground-state Ag and Sg baryons, one 
defines two hadronic parameters, Ai and A2, as 

{B{v) I K I B{v)) = Ai, 

{B{v) \ S \ B{v)) = dM X2. (2) 

where (Im is zero for Aq and — |, 1 for Sg, Sg baryons, respectively. The constant (Im 
characterizes the spin-orbit interaction of the heavy quark and the gluon field. Therefore, 
the mass of heavy baryon up to order l/mq corrections can be written in a compact form 

M = mQ + A- -^(Ai + dAfA2) , (3) 

where parameter A is the energy of the light degrees of freedom in the infinite mass limit. 
Thus the splitting of the spin 1/2 and 3/2 doublets is 

Sq^-S^ = ^A2. (4) 

The hadronic parameters Ai and A2 are nonperturbative ones that should be either 
determined phenomenologically from experimental data or estimated in some nonpertur- 
bative approaches. A viable approach is the QCD sum rules P] formulated in the framework 
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of HQET |]7[. This method allows us to relate hadronic observables to QCD parameters 
via the operator product expansion (OPE) of the correlator. In the case of heavy mesons, 
those two matrix elements thus masses had been calculated to perfect first by Ball and 
Braun and latterly by Neubert P] taking a different approach. Masses of excited meson 
states had been calculated up to l/mq order in [jlO[. For the case of heavy baryons, there 
are several attempts to calculate the baryonic matrix elements of /C and S. Using HQET 
sum rules, Colangelo et. al. have derived the value of Ai for Aq baryon |[Tl|]. Furthermore, 
the baryonic parameters Ai and A2 for the ground state baryons had been calculated in |]T2| 
by evaluating the two-point correlation functions. The mass parameters of the lowest lying 



excited heavy baryons had also been determined recently in [r^|. In the present work we 
shall calculate the baryonic parameters Ai and A2 for ground state Ag and Sg baryons using 
QCD sum rules in the HQET. Following Ball et. al. [§ and Neubert's |^ work done for the 
meson, we adopt these two approaches, named as direct approach and covariant approach, 
to evaluate the three-point correlators and obtain the values of baryonic parameters. It is 
of interest to compare the two methods in the analysis. 

The remainder of this paper is organized as follows. In Sec. [11 A| we introduce the 
interpolating currents for baryons and briefly present the two-point sum rules. The direct 



Laplace sum rules analysis for the matrix elements is presented in Sec. |II B| . Another 
feasible approach (covariant approach) to this aim can be found in Sec. 11 C Sec. |IT1| is 
devoted to numerical results and our conclusions. Some comments are also available in 

Sec. m 



II. DERIVATION OF THE SUM RULES FOR Ai AND A2 

A. Heavy baryonic currents and two-point sum rules 

The basic points in the application of QCD sum rules to problems involving heavy 
baryon are to choose a suitable interpolating current in terms of quark flelds and to deflne 
the corresponding vacuum-to-baryon matrix element. As is well known, the form of inter- 



polating currents for baryon with given spin and parity is not unique plHlq] , the choice of 
which one is just a question of predisposition. The most generally used form of the heavy 
baryon current can be written as 



r = eaUqI''CTTql)T'hl, (5) 

in which C is the charge conjugation matrix, r is the flavor matrix which is antisymmetric 
for Aq baryon and symmetric for Eg baryon, F and F are some gamma matrices, and a, 
b, c denote the color indices. F and F' can be chosen co-variantly as 
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for Aq baryon, and 



for Eq baryon, and 



75 , r' = 1 , (6) 



7m , r' = (7^ + v^) 75 , (7) 



11 2 

r = 7^, T' = -gf,^ + -'y^'y^--{'yf,v^-'j^v^) + -v^Vf,, (8) 

for Sq baryon. Also the choice of F is not unique. We can insert a factor before F defined 
by equations The currents given by Eqs.®-® are denoted as and that with ^ 

insertion as jg, which are two independent current representations. 
The baryonic couphng constants in HQET are defined as follows 

{0 \ f \ A{v)) = Fj,u, 
(0 I f I = Fj,u, 

(0|/|S*(^;)) = i=Fs.u^ (9) 

where u is the spinor and Ua is the Rarita-Schwinger spinor in the HQET, respectively. 
The coupling constants Fs and are equivalent since Eg and Eg belong to the doublet 
with the same spin-parity of the light degrees of freedom. 

The QCD sum rule determination of these coupling constants can be done by analyzing 
the two-point function 

t J dxe''-'{0 I T{f{x)r{0)} I 0) = l±lTr[rr+]U{u), (10) 

where k is the residual momentum and uj = 2v ■ k. It is straightforward to obtain the 
two-point sum rule: 

p2 ^-2A,/T ^ 911 SrM/T) - {'^G') SMJT) + i imf. (11) 
The functions 6n{uJc/T) arise from the continuum subtraction and are given by 

(5„(x) = - / dtre-* = 1 - e-^' ^ |- . (12) 

The second term of the last equation is assigned to the continuum mode, which can be 
much larger than the ground state contributions for the typical value of parameter T due 
to the high dimensions of the spectral densities. 
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B. The direct approach 



In order to evaluate the matrix elements Ai and A2 we consider the three point corre- 
lation functions with K, and S inserted directly between two interpolating currents at zero 
recoil as below 

^' jdxj dye''-^-'''-y{0 \ T{j^{x)ICiO)]^iy)} \ 0) = ^Tr[TT+] Tk{uj,uj') , 

^^ jdxj dye''-^-'''-y{0 \ T{j^{x)S{0)]^ (y)} \ 0) = dM^Tr[Tr^] Ts{uj,uj') , (13) 

where the coefficients Tk{uj,uj') and Ts{uj,uj') are analytic functions in the "off-shell ener- 
gies" u = 2v ■ k and u' = 2v ■ k' with discontinuities for positive values of these variables. 
Saturating the three-point functions with complete set of baryon states, one can isolate 
the part of interest, the contribution of the lowest-lying baryon states associated with 
the heavy-light currents, as one having poles in both the variables u and u' at the value 
uj = uj' = 2A: 



Tk{uj,uj') 



(2A-cj)(2A-cj' 



where the ellipses denote the contribution of higher resonances. In the theoretical calcula- 
tion of the correlator it is convenient to choose the residual momenta k and k' parallel to 
the V, such that A;^ = ^v^ and k'^ = ^v^. 

The leading contribution to the matrix element of kinetic energy is of order 1, whereas 
to the chromo-magnetic interaction is of order a^. Confining us to take into account these 
leading contributions of perturbation and the operators with dimension D < 6 in OPE, 
the relevant diagrams in our calculation are shown in Fig. 1. and Fig. 2. The calculation 
of the diagram (a) in Fig. 2. is the most tedious one. It can be computed using Feynman 
parameterization and the integral representation of the propagators, which is the standard 
technique |18|,|19[. The factorization approximation has been used to reduce the four-quark 
condensates to (gg)^ in the calculation. 

On theoretical side the correlators Tk{oj,uj') and Ts{uj,uj') can be casted into the form 
of integrals of the double spectral densities as 

T^(o;,o;') = //^^Px(s,s'), 

T.(^,-') = //j^^P5(s,s'), (15) 
where the double spectral density functions are 
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TC 



p|(s,s/) = -^[e(s-s') r dx(s-x)(s'-x)a;3 + e(s'-s) dx{s - x){s' - x)x^] 

Z'^TT^ Jo Jo 

^ :{^G') s'6{s - s') + ^{qqf [s6{s') + s'S{s)]. (16) 



487r2 ' TT ' ' 'Svr 
The unitary normalization of flavor matrix Tr[TT~^] = 1 has been applied to get those 
densities. Here we use the numbers 1, 2 to denote the results corresponding to the different 
choice of currents and jg. Those we do not discriminate with numeric superscripts 
meas that with or without insertion the results are identical. Following Refs. PD|- P^] , 
we then introduce new variables = ^{uj + u') and uj_ = oj — cu', perform the integral 
over and employ quark-hadron duality to equate the remaining integral over uj+ up 
to a "continuum threshold" Uc to the Borel transform of the double-pole contribution in 
(|l^. Then following the standard procedure we resort to the Borel transformation i?^, 
B^, to suppress the contributions of the excited states. Considered the symmetries of the 
correlation functions it is natural to set the parameters r, t' to be the same and equal to 
2T, where T is the Borel parameter of the two-point functions. We end up with the set of 
sum rules 

-4Af ■ V^e-^'/- = - ( *'''^'/^>' 
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a. 



4A,F2e-2^-/^ = l^^T«57(^c/T) - ^(^G^) 5,{ujjT) + ^^^(gg)^ 5i(..,/T), (17) 
vr^ vr ovr^ vr Svr 

It is worth noting that the next-to-leading order as corrections have not been included in 

the sum rule calculations. However, the baryonic parameter obtained from the QCD sum 

rules actually is a ratio of the three-point correlator to the two-point correlator results. 

While both of these correlators are subject to large perturbative QCD corrections, it is 

expected that their ratio is not much affected by these corrections because of cancellation. 

On the other hand, we have only calculated the diagonal sum rules by using the same type 



5 



interpolating current in the correlator. As to the non-diagonal sum rules, the only non- 
vanishing contributions in the OPE of correlator are terms with odd number of dimensions, 
thus the perturbative term gives no contribution. The resulting sum rules are dominated 
by the quark-gluon condensates. It is expected that the non-diagonal sum rules will give 
no more information than diagonal ones. This has been proved to be true in the analysis 



of Ref. \ll 



C. The covariant approach 

In the previous subsection we have completed the task of the determination of the 
matrix elements for both the operators of kinetic energy and chromo-magnetic interaction 
by direct calculation of three-point correlation functions. In fact, there exists a field- 
theory analog of the virial theorem p3|,p^ in consideration of the restrictions the equation 



of motion and the heavy quark symmetry imposing on baryons, which relates the kinetic 
energy and chromo- interaction to each other and ensure the intrinsic smallness of the kinetic 
energy explicitly. In this subsection we shall follow Neubert's procedure and take those 
restrictions into account to deduce a new result of the kinetic energy (the chromo-magnetic 
interaction is identical). 

The main idea of that procedure is that the coefficients of the covariant decomposition of 
the bilinear matrix elements, the so called invariant functions, can be related to the kinetic 
energy and chromo-magnetic interaction at the zero recoil. Following the discussion in 
P,p5| we have the general decomposition (see Appendix): 



(A I Ka^.ig.G^^'K, \ A') = 0i {v\ v) {v'^v, - v^v',)ua^''u\ (18) 
for Aq baryon, in which u is the spinor in HQET, and 

(E I Ka^^igsG^^^K \ S') = 0^^(t;', v)^!''a^^m'f' , (19) 
for Sq baryon, where bears the decomposition 

+ (piiv'^Vy - v'^v^,)ga(5 + hW^Vy - vlv^)v'^vp, (20) 

in which we use the covariant representation of the doublets \E'^ = + ■:^{Vfj, + 7^)u with 
restrictions -^u = u, v^u^ = and 7^^^^ = 0. The normalization of those coefficients at 
zero recoil is 0i(l) = — |Ai for Aq baryon and 

±A2 = 20i(l), 

±Ai = 0o(l) = 0i(l) - 2(02(1) - 03(1)) - 304(1), (21) 
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for Tjq baryon. The foregoing minus sign corresponds to the Sq baryon and plus to Tiq 
baryon. 

Let us now derive the Laplace sum rules for the invariant functions (j)i{w). The analysis 
proceeds in complete analogy to that of the Isgur-Wise function. We shall only briefly 



sketch the general procedure and refer for details to Refs. [20,21|. We consider, in the 



HQET, the three-point correlation function of the local operator appearing in ([19|) with 
two interpolating currents for the ground-state heavy baryons: 

^2 j dxdi/e^'=-^-*^'-^(0|T{/(x),M(7,rG"^X.(0),?' (i/)}|0) 

= Kp{v',v,k',k)T'^^-^T^-^T',^ (22) 

where k and k' are the residual momenta. The Dirac structure of the correlation function, 
as shown in the second line, is a consequence of the Feynman rules of the HQET. <l>(^^ obeys 
a decomposition analogous to (pO|), with coefficient functions cj', w) that are analytic 
in the "residual energy" uj = 2v ■ k and u' = 2v' ■ k', with discontinuities for positive values 
of these variables. These functions also depend on the velocity transfer w = v ■ v' . 

The lowest-lying states are the ground-state baryons B{v) and B'{v') associated with 
the heavy-light currents. They lead to a double pole located at u = u' = 2A. The residue 
of this double pole is proportional to the invariant functions (piiw). We find 

^"'"'"^= (c.-2A)(..--2A) ' ^22^ 

where Sc is the structure constant, 1 for Aq, — for Sg and 1 for Eq baryons. In the 
deep Euclidean region the correlation function can be calculated perturbatively because of 
asymptotic freedom. Following the standard procedure, we write the theoretical expressions 
for $j as double dispersion integrals and perform a Borel transformation in the variables u 
and u', then set the associated Borel parameters equal: t = r' = 2T. All goes like that in 
the direct approach, we introduce new variables = ^{u + oj') and U- = u — uj' , perform 
the integral over tu., and get the Laplace sum rules at zero recoil: 

8s^F'<p,e-'^-/^ = --^^T%{ujjT) (24) 



for Eq baryon, and 



SFVie-^^^/^ = --^^T%{u;jT), (25) 



TT'^ TT 
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for Aq baryon. After some simple algebra we find 
for Sq baryon, and 

- 4AiF^e-2^^/^ = -^^<5,(.;,/T), (27) 

for Aq baryon. The minus sign of the Ag baryon result may seem bizarre, in Sec.|T^ we 
will return to dwell on this point. 



III. NUMERICAL RESULTS AND CONCLUSIONS 



In order to get the numerical results, we divide our three-point sum rules by two-point 
functions to obtain Ai and A2 as functions of the continuum threshold ujc and the Borel 
parameter T. This procedure can eliminate the systematic uncertainties and cancel the 
parameter A. As for the condensates, we adopt the standard values 

(gg) = -(0.23 GeV)^ 

rv 

{^G^) = 0.012 GeV^ (28) 
vr 

From two-point sum rules one has known that there exist stable windows between 0.8 < 
T < 1.2 GeV and = (2.2 - 2.7) GeV for Aq baryon, and 0.7 < T < 1.1 GeV and 
cjc = (2.6 — 3.3) GeV for Sg baryon. The stability window for three-point functions starts 
almost from values of the Borel parameter at 0.7 GeV and stretches practically to T 00. 
It is known that stability at lager Borel parameter could not give any valuable information, 
since in this region the sum rule is strongly effected by the continuum model. The usual 
criterium that both the higher-order power corrections and the continuum contribution 
should not be very large restricts the working region considerably. In the case of the three- 
point functions, the results are severely smeared by the continuum contributions for the 
high dimension of spectral densities and thus it is very difficult to ensure the contribution 
of the continuum mode is small. The working region of the three-point functions should 
be determined by the stable region of the two-point functions. So it does not necessarily 
coincide with the stable windows of the three-point functions 0|. Thus we find our working 
region for three-point functions is T = 0.8 — 1.2 GeV for Aq baryon and T = 0.7— 1.1 GeV 
for Sq baryon. The results for Ai of the direct approach with two different choices of 
currents are shown in Fig. 3 for Aq baryon and Fig. 4 for Sq baryon. Results for kinetic 
energy of Aq and Sg baryons obtained by the covariant approach are presented in Fig. 5. 
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The forms of the chromo-magnetic interaction obtained by both approaches do not differ 
from each other, so we plot that unique curve in one figure, Fig. 6. For the Ag baryon we 
obtain the residual mass Aa = 0.8 ±0.1 GeV and 

-Xl = 0.4 ±0.1 GeV^ 

-A? = 0.5±0.1 GeV^, (29) 

in the direct approach, where the superscripts denote the different choices of currents, and 

- Ai = -(0.08 ± 0.02) GeV^ (30) 

in the covariant approach. For the Eg baryon the effective mass we obtained is As = 
1.0 ± 0.1 GeV and 

-X\ = 0.7 ±0.2 GeV^ 

-A^ = 1.0 ± 0.2 GeV^ (31) 
in the direct approach, where the superscripts also denote different currents, and 

- Ai = 0.11 ±0.03 GeV^ (32) 
in covariant approach. For the chromo-magnetic interaction for S baryon the results read 

A2 = 0.23 ± 0.02 GeVl (33) 
Then we get the splitting of the spin 1/2 and spin 3/2 doublets is 

- = ^ A2 = 0.35 ± 0.03 GeV^ (34) 

All error quoted before is due to the variation of Borel parameter T and the continuum 
threshold Uc- When it is scaled up to the bottom quark mass scale there will be a factor 
~ 0.8 approximately due to the renormalization group improvement. 

As for the effects on the correlation function of the different choices of the interpolating 
currents we may assert some facts and inspections. From the preceding numerical results 
it is clear that the interpolating currents with the |^ insertion give a considerable larger 
result to the kinetic energy than those without the insertion. Nevertheless, the two-point 



sum rules do not differ with the different currents [|T^,|T^. From our calculation it is 
explicit that the sum rules associated with chromo-interaction insertion are identical. In 
our covariant calculation we find that the invariant functions do differ from each other 
generally, but interestingly they coincide at the zero recoil thus the chromo-interaction 
and kinetic energy do not take a different form. Naively, we can tell that the disparity of 
the two forms of the kinetic energy obtained in our direct calculation mainly comes from 
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the Lorenz structural differences of tlie two interpolating currents. It may be noted that 
derivative operator acts differently on the currents with or without |^ insertion, thus with the 
insertion it is easier for the continuum contamination to go into the correlation functions. 
It is urgently needed to exclude the continuum contribution which smeared heavily the 
results of both the direct and covariant approaches. It is this continuum contamination 
that makes the prediction of the kinetic energy more intriguing. All previous theoretical 
calculations with QCD sum rule approach or lattice calculation give various results and 
can differ from each other by several times p!l|JT^ J^,P,p^,p^ . Current experimental data is 



not enough to judge which one is right and what we can get is some restrictions on the 
kinetic energy |^ or a rough estimate of the kinetic energy extracted from experimental 



data with some assumption [Q. As demonstrated in Refs. p^j29[] using a toy model of 
harmonic oscillator, the main origin of the discrepancy between the direct and covariant 
approach is the continuum smeared contribution. In the direct approach the first excited 
contribution plays an important role. If we want to suppress this contribution, we go to 
such a large Borel parameter that the power corrections blow up. For acceptable Borel 
parameter, we get an over-estimated sum rule for the the kinetic energy. In the covariant 
approach (via virial theorem) the situation is especially bad. The excited contribution 
consists of two components - the diagonal transitions and off-diagonal ones - and each one 
is large, but they have opposite relative sings. For highly excited states the sign-alternating 
terms are smeared to zero after summation. However, the first two terms do not cancel 
with each other and screen the ground state contribution. Thus a lower-estimated result 
will be obtained. The minus sign before the kinetic energy of Ag baryon in the covariant 
approach may be seen as a manifestation of this assertion. Due to the unknown weight, we 
cannot annihilate those contributions by weighted averaging just like that in the Quantum 
Mechanics. But we may safely take the results of the direct and covariant approach as 
lower-bound and higher-bound of the kinetic energy parameter, respectively. Then, follow 
pOj , to take the mean value of the direct and covariant approach results as an rough 



estimate. The result thus obtained is 



A} ~ 0.18 ±0.06 GeV^ 

-A^ ~ 0.24 ±0.06 GeV^ (35) 



for Aq baryon and 



- Ai ~ 0.39 ±0.12 GeV^ 

-A? ~ 0.54±0.12 GeV^ (36) 

for Eq baryon. Taking all results obtained the mass of the ground state baryon is on hand. 
From and ttia^ pT|, we determine the heavy quark masses rric — 1.41 ± 0.16 GeV 
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and nib — 4.77 ± 0.12 GeV. In the determination we have taken the average of the results 
obtained from two interpolating currents to be the physical pole masses of the heavy quarks 
because that the difference of the corresponding mass does not exceed the error bar. These 
values give the following results: 



with interpolating current jl and 



ms, ^ 2 


47 ±0.20 GeV, 


ms. ^ 2 


59 ± 0.20 GeV, 


msfc ^ 5 


79 ±0.13 GeV, 


m^* — 5 


82 ±0.13 GeV, 


and 




~ 2 


52 ± 0.20 GeV, 


ms* ^ 2 


64 ± 0.20 GeV, 


ms^ ~ 5 


80 ±0.13 GeV, 


m^* ~ 5 


83 ±0.13 GeV, 



(37) 



(38) 

with interpolating current j^. The spin average of the doublets are free of the chromo- 
interaction contribution and thus free of the uncertainties involved in the calculation of A2. 
Average over the doublets we have the quantity 

^ (Ms^ ± 2Ms^) = mg ± As ± 

which is more reliable. For the c quark case, it is 2.55 ± 0.20 GeV with current jl and 
2.60 ± 0.20 GeV with current j^. For the b quark case it is 5.81 ± 0.13 GeV with current jl 
and 5.83 ± 0.13 GeV with current £. Experimentally Ms, = 2453 ± 0.2 MeV There 



is experimental evidence for S* at Ms* = 2519 ± 2 MeV |^2[. If we take this value for 
E*, we have |(Ms^ ± 2 M^*) = 2497 ±1.4 MeV which is in reasonable agreement with 
the theoretical prediction. For lack of experimental data the corresponding quantity for 
the bottom quark will be checked in the future. If we take the preceding masses of the 
charmed S baryons the splitting thus reduced is 0.33 GeV^ and our theoretical splitting is 
in considerable agreement with the experimental data. 

As the kinetic energy of Aq baryon can be related to the spectrum via the kinetic 
energy of meson as p3| n 

(mA, - mAj - {mB - niD) = [Ai(A,) - \,{B)] - ^) ± 0(l/m|) , (39) 



^The relation between /i^ in Ref. |33] and Ai in this paper is //^ = — Ai 
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where = i(m^ + 3m^) and tTid = j{m£) + 3m^) denote the spin- averaged meson 
masses, the difference between the Icinetic energy of B meson and that of A;, baryon can 
be extracted as 

Ai(Ab) - Xi{B) = 0.01 ± 0.02 GeV^ , (40) 



which is consistent with the value obtained in Ref. [^]. Resorting to the recent experimen- 
tal data for the mesonic kinetic energy parameter obtained in the inclusive semileptonic B 
decays [Q, — Ai = 0.24 ±0.11 GeV^, one can thus get the value of baryonic kinetic energy 



as 

- Ai(A;,) = 0.23 ±0.13 GeV% (41) 

which is in reasonable agreement with our theoretical prediction given in (^). 

For conclusions, we have calculated the l/mg corrections to the heavy baryon masses 
from the QCD sum rules within the framework of the HQET. Two approaches have been 
adopted in the evaluation of the three-point correlators. Our final results read 

= mg ± As ± 77^(0.16 ± 0.12 GeV^), 

Ms* =mQ± As ±77^(0.51 ±0.12 GeV'), (42) 

for interpolating current without insertion and 

Msq = mQ ± As ± 77^(0.31 ± 0.12 GeV^), 

Ms* =mQ±As±-^(0.66±0.12GeV2), (43) 
« 2mQ 

for interpolating current with insertion. The l/mg corrections are small. We have 
taken the mean value of the direct and covariant approach as the rough estimate of the 
kinetic energy parameter Ai. Our theoretical predictions are in agreement with the recent 
experimental data. For a more precise treatment of the kinetic energy, more sophisticated 
technique to distinguish the smearing continuum contribution is in urgent necessity to be 
developed. 
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APPENDIX A: THE DECOMPOSITION OF THE BILINEAR MATRIX 

ELEMENT 

In this appendix, we present the decomposition of bihnear matrix element. There exists 
the decomposition of Aq baryon, we present it here merely for completeness and convention. 
First, let us consider the bilinear matrix element over Aq baryons 

(A I K{-i%)T^'iD,K, I A') = ^^,{v\ v)uT^''u\ (Al) 

the coefficients obey the symmetric relation ip^i,{v', v) = ipl^iv, v'). It is convenient to write 
the coefficient into the sum of symmetric and anti-symmetric parts tf^^i, = \ [ip-^j, + ip^^] 
which can be presented covariantly as 

= ^hnu + i^liv + v')^,{v + v% + i^i{v - v')^{v - (A2) 

the HQET equation of motion implies that v'^^ip^i, = from which we can obtain the 
relations between those coefficients 

V'f + (1 + y)i^l + (1 - y)i^l + yi^t = 0, 

+ + (2/-l)^f-V'i' = 0, (A3) 

with 

hi%Vh' + hiD^Vh' = id^{hVh') (A4) 

bear in mind we can get 

(A I Ji^mD^iD.h^. I A') = ip^,{v', v)ur^''u' + A{v' - v)^^,ur''''u' , (A5) 

usmg X dependence of state in HQET | B{x)) = e"^^^'^ | 5(0)). The are defined as 

(A I KT'^iD^K, I A') = ^^^uV^u', (A6) 

Similarly, we can define the matrix elements for the operators of kinetic energy and chromo- 
interaction over baryon states with different velocities 

(A I Ka^JC^^K^ I A') = Mv'^^v, - v^vDua'^'^u', (A7) 

(A I K{iD^fVK, I A') = <pouVu\ (A8) 
once such defined, the ■^j can be expressed via two 0^ 
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y + 1 

^1 = 00 + #1 + A'e 



y + r 



2(y-l) 2(y + l) ^' 

= (A9) 

1 + 1/ 

the normalization of 0o, 0i are 00(1) = Ai, 0i(l) = — |0o(l), thus we get the desired resuh. 

Generahzation can be made to the higher spin states such as Sq baryon. The procedure 
goes almost the same. The only difference lies on the decomposition of the matrix element. 
Hence we will give the forms of decomposition and the final desired relation, others we will 
not dwell on. The covariant representation of the doublet is ^'^ = + -^{v^ + 7^)m. The 
matrix element is 

(E I K(-i^^)mD,K, I E') = ^'^^{v', v)^^r^''%, (AlO) 

in which the coefficients obey symmetric relation -0^^ = ip^^{v',v). Adopt the same 
symmetric and antisymmetric decomposition of the coefficients like that in Aq baryon case, 
we have 



+ i^Udat^v^vp + g,,av^vp + g^^v'^v'^ + g^pv'y^) + ipl {v' - v)f,{v' - v)WaV^ 

+ 0'6 - - v)^gap + i^r {v' + v)i^{v' + v)^v'^vp + t/jl {v' + v)i^{v' + v)^gai3 

+ ^I^Mi^v, (All) 



(S I Kr{-iD,)K> I E') = lt^{v\ v)^j:^^%, (A12) 



introduce other universal parameters in the leading order 

as usual, i'^^{v,v') can be decomposed into the general form 

C^(^^, v' = ii{v + v')^ga/3 + 6(^' - v)^gap + ^siv + v')^vpv'^ 

+ ^4{v' - v)^Vf}v'^ + ^5v'a9(^u + ^eVf^g^i., (A13) 

the equation of motion implies that v"^^^^ — and v"^i/j^^ — from which we can derive 
relations 
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5 



wi)^ + v^f + (1 - w)^f + (1 + w)i), 



i)^ + «7V^5^ + + (1 - t/7)^5^ + (1 + + i)l = 0, 

+ {w- - - + (1 + = 0, 



and 



take the difference of the two terms in Eqs. ( [A12| ) and use ( |A4| ) we can reach 



(A14) 



(i + ^)6 + (i-w^)6 = o, 

{l+w)^s + {l-w)U + ^e = 0, (A15) 



6 = ^^A, 

_ W - 1 C2 . 

w + 1 2 

fs = &, (A16) 

where Ci , C2 parameterize the matrix element 

(S I /i,r/i„, I S') = (ci(7«^ + C2t;/30^"r^'^. (A17) 

The matrix elements for the kinetic energy and chromo-magnetic interaction are defined 
similar to those for the Ag baryon 

(S I KitD^fTK, I S') = (0o^7a/3 + 00^^/30 ^'"r^'^ (A18) 

(S I /i.a^.^G^^/i,, I S') = ^^^^^a^'^^lj, (A19) 

where bear the same decomposition as ip"^^ and they have simple relations between 
each other 

01 = i^t, 

02 = ^2 - ^eA, 

03 = V^3''-e6A, 
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04 = il^t - S^iA, 

05 = i't - 2e3A, 

00 = 2^1 + + (1 - w)r^ + (1 + + 2(1 - w)i2K 

00 = SV'I + (1 - w)V'5 + (1 + w)V'7 + 2(1 - w)iiK, (A20) 

the normalization condition is that 0i(l) = 0o(l) = B\i where A is —1/2,1/2 and 
S is 1, —1 for Eg, Eq respectively. At zero recoil Ai can be expressed via 0o 

0o(l) = 0i(l) - 2[02(1) - 03(1)] - 304(1) . (A21) 
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Figure Captions 

Fig. 1. Non- vanishing diagrams for the kinetic energy : (a) per- 
turbative contribution, (b) to (e) gluon-condensate. The kinetic en- 
ergy operator is denoted by a white square, the interpolating baryon 
currents by black circles. Heavy-quark propagators are drawn as 
double lines. Diagrams (b) to (e) are calculated in Fock-Schwinger 
gauge. The lower right vertices of those diagrams are set to the 
origin in coordinate space. 

Fig. 2. Non-vanishing diagrams for the chromo-magnetic inter- 
action: (a) perturbative contribution, (b) gluon-condensate, (c) 
quark-condensate. The chromo-magnetic interaction (velocity- 
changing current) operator is denoted by a white square, the in- 
terpolating baryon currents by black circles. 

Fig. 3. Sum rules for Aq baryon: (a) for (b) for j|. The 
dash-dotted, dashed and solid curves correspond to the threshold 
ujc — 2.4, 2.6, 2.8 GeV, respectively. The working region is T = 
0.8- 1.2GeV. 

Fig. 4. Sum rules of the kinetic energy for Eq baryon: (a) for jl, 
(b) for jl. The dash-dotted, dashed and solid curves correspond 
to the threshold Uc — 2.9,3.1,3.3 GeV, respectively. The working 
region is T = 0.7 - 1.1 GeV. 

Fig. 5. Covariant sum rules of the kinetic energy: (a) for Aq 
baryon, (b) for Sq baryon. The dash-dotted, dashed and solid 
curves correspond to uJc = 2.4, 2.6, 2.8 GeV for Aq baryon, and 
ujc = 2.9,3.1,3.3 GeV for Eg baryon. The working region is T = 
0.8 - 1.2 GeV for Aq baryon and T = 0.7 - 1.1 GeV for Eq baryon. 

Fig. 6. Sum rules for the chromo-magnetic interaction. The 
dash-dotted, dashed and solid curves correspond to lOc — 
2.9, 3.1, 3.3 GeV. The working region is T = 0.7 - 1.1 GeV. 
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FIGURES 
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